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ABSTRACT 

Motivated by recent Belle measurement of the branching ratios — »• 
TT+TT^i^" and B~ — ^ K^K~K~ and B^ XcqK^., we investigate the CP 
violating asymmetry in the partial widths for the decays B^ -k^-k^ K~ 
and B^ K^K^K~ , which results from the interference of the nonresonant 
decay amplitude with the resonant amplitude B^ XcoK^ n^n^K^ and 
B^ — > XcqK^ — ^ K^K~K~ . By taking 7 ^ 58° we predict that the partial 
widths asymmetry is reaching 10% for the B^ —>■ tt'^'k~K~ decay and 16% 
for the B- K+R-R- decay. 



The extraction of the CP violation phases a, (3 and 7 has stimulated 
many studies during last decade ((see e.g. [|l|). Among them the phase 7 
is most difficult to constrain. Numerous proposals were aimed to determine 
its size (see e.g. [0]). The partial width asymmetry in charmless three-body 



decays offers a chance to get more information on sin 7 p[| - 110 . 

The basic idea was that the nonzero asymmetry results from the inter- 
ference of the nonresonant three-body decay amplitude and the resonant 

— > XcoTi"^ 7r+7r^7r~ decay amphtude 0, ||, |10[- AH these approaches 
were based on a theoretical prediction [§] for the decay rate of B~ Xcot<'~ 
which is difficult to obtain in the factorization model. Instead of relying on 
this prediction for this decay rate, we can now use the recent Belle collabo- 
ration measurement of the decay rate |TT|, IT2I 



BR{B- ^ XcoK-) = (Q.Otii) X 10-^ (1) 

which is surprisingly large, comparable to the B^ J/ipK" decay rate. In 
addition, the Belle collaboration has anounced the measurement of the decay 
rates 

BR{B- K+K-R-) = (37.0 ± 3.9 ± 4.4) x 10"^ (2) 

and 

BR{B- n+n-R-) = (58.5 ± 7.1 ± 8.8) x 10"^ (3) 

Although these measurements do not distinguish yet the contributions to 
the branching ratios of the nonresonant decays, they constrain their mag- 
nitudes. On the other hand, the charmless two-body decays of B mesons 
are stimulating many research projects. The factorization approximation 
is usually used and recently great improvement has been made in the un- 
derstanding of -B ^ TTTT and B —>■ ttK decay modes (e.g. ||1^). The de- 
cay mechanism of the nonleptonic charmless three-body B decays is much 
less understoood. The intermediate resonance states might be simpler to 
understand due to better knowledge of the two-body B decays. Due to 
the lack of understanding of the three-body nonresonant decay modes of 
B mesons, we will adopt the existing approaches The first is the 

use of the factorization approach. Then we reduce the problem to the cal- 
culation of the matrix elements of currents (or density operators). In the 
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evaluation of these matrix elements, the main contribution to the nonres- 
onant MMK^, M = i^+,7r+, amplitude comes from the prod- 

uct < MM\{ub)v~A\B^ > < K^\{su)v-a\0 > where {qiq2)v-A denotes 
gi7M(l - 75)g2, or from Eg < MM|(g(l - Ts)^!^- > < K-|(s(l + 75)g)|0 >. 
For the calculation of the matrix element < MM\{ub)v-A\B~ > we extend 
the results obtained in where the nonresonant K~tt^Iv decay was 

analyzed. In this analysis the experimental result for the branching ratio of 
the nonresonant K^tt^Iu decay was successfully reproduced within a 

hybrid framework which combines the heavy quark effective theory (HQET) 
and the chiral Lagrangian (CHPT) approach. 

The combination of heavy quark symmetry and chiral symmetry has been 



quite successful in the analysis of D meson semileptonic decays []T3[ - [pi|]. 
Heavy quark symmetry is expected to be even better for the heavier B mesons 
T^. However, CHPT could be worse in B decays due to the large energies 



of light mesons in the final state [|5[. In order to adjust the model to be 
applicable in the whole kinematic region we retain the usual HQET Feynman 
rules for the vertices near and outside the zero- recoil region, as in |T8|, p!9|] , 
but we include the complete propagators instead of using the usual HQET 



propagator |jT^, pT], p2 



In the following we systematically use this model [Tj, ^ to calculate the 
nonresonant — >■ MMK^, M = K^,7i^, decay amplitude. We find that 
contrary to the previous cases [§] penguin contributions are very important. 

First we discuss the contributions to the nonresonant decay amlitude and 
then we discuss the partial width asymmetry. 

The effective weak Hamiltonian for the nonleptonic Cabibbo-suppressed 
B meson decays is given by 



Heff = ^[K;K6(CiOi„ + C2O2J + V:MciOi, + C2O2C)] 
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- Y.i[yubv:,cr + v,bV:A + Vt,v;A)o^] + h.c. (4) 

1=3 

where the superscripts u, c, t denote the internal quark. The tree-level op- 
erators are defined as 

Oi, = g7M(l-75)&^7''(l-75)g, (5) 

with q = u,c and 

02g = g7/.(l- 75)9^7^1 -75)6. (6) 
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We rewrite O3 — Oq, using the Fierz transformations, as follows: 



03= E s-f,il-^5)bqr{l 



q=u,d,s,c,b 



O4 



J2 s7/.(i -75)gg7''(i -75)^, 

q=u,d,s,c,b 

^7m(i -75)&g7''(i + 75)g, 

q=u,d,s,c,b 

-2 J2 s(l + 75)gg(l-75)6, 

q=u,d,s,c,b 



(7) 

(8) 
(9) 
(10) 



The operators 07,...,0io denote the electroweak penguin operators. Due to 
the smallness of their Wilson coefficients we do not include them in our 
analysis [H ||]. 

The factorization approximation is obtained by neglecting in the La- 
grangian terms which are the product of two color-octet operators after Fierz 
reordering of the quark fields. The effective Lagrangian for non-leptonic de- 
cays are then given by @ with Cj replaced by a^. In our calculations we 
use next-to-leading Wilson coefficents (e.g. O)- In @) 0; ^ it was 
found that 03 and 05 are one order of magnitude smaller than and ag and 
therefore we can safely neglect the contributions from O3 and O5 operators. 
3, 



For A^, 



5 GeV, we use [26 



ai 



1.05, 

-0.043 -0.016Z, 
-0.054 -0.016Z. 



a2 



0.07, 



(11) 



It is important to notice that the imaginary parts of the Wilson coefficients 
for the penguin operators are due to the internal charm- and up-quark loop 
exchanges. They introduce additional strong phases, affecting the CP asym- 
metries. 

The use of factorization gives for the matrix element of Oi operator for 
the nonresonant tt^ti^K^ decay, following the results of H 



<7r+(pi)7r {p2)K{p3)\Oi\B (pb) >nr 
< K~{p3)\{su)v^a\0 >< Tr'{p2)n+ipi)\{ub)v-A\B'{pB) >nr ■ 



(12) 
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To evaluate the matrix element < 7i~ {pi)7i^ {p2)\{ub)y^A\B^ (ps) >nr we 
will also use the results obtained previously in the analysis of the nonres- 
onant -k^K'Ivi decay width [Q. We write the matrix element 

< Ti^ {pi)T^^{p2)\{uh)v-A\B~ {jpb) > in the general form 

< 71' {pi)7i^ {p2)\iub)v-A\B' (pb) > = ir{pB - P2 -Pi)fi 

+iw+{p2 + Pi)^ + iw.{p2 - Pi)^, - 2he^a(3jPBP2Pi ■ (13) 



The form factors and r"*" for the nonresonant decay are given in 0, ^, |T4 
together with a detailed description of our hybrid model. Using the following 
set of Mandelstam's variables: s = {pb — PsY, t = {pb — pi)^ and u = 
{Pb — P2)'^ the formf actors can be written as 



wT{s,t) 



r 3/2 1/2 / 2 2 j-\ 



2 

B* 



/1/2 {t - 2m 
+ 7Tl^-#f^A(2t + .-m|-m^-2mD, (14) 



2/1/2 /1/2 2m 



w'_ 



{s,t) 



r 3/2 1/2 

9 fB*mj^^mj^ 



+ 



/1/2 t 



m 



■[1 + 



B* 



2mL 



(15) 



r"^(s,t) 



X 

+ 
+ 



1+131 



{2t + s 



nil 



/1/2 2 
^B WfB"m*j^imB 



m 



2ml] 



ImB' 



B' 



m B ml — m\i 



mB" 



/1/2 



mn — m 



B" 



Us -ml- ml) - j-^{t + m\- ml){m\ -ml- t)\ 



t — m 



B' 
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'mB f. 



B'- 



/1/2 V "^B' ^ 



[m 



m 



B 



mi 



t) 



B'* 



B 



+ 



m^a2 i 
/1/2 2m| 



(2t + s — m\ — m\ — 2ml] 



(16) 
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With the use of these resuhs we write down 

< K-{p3)'K+{pi)7T-{p2)\Oi\B{pb) >nr = 



i/3(s + 2t - m| - 2ml ' ^D^-T (17) 



The parameters ai^2 are defined in pT[, g is the B*Btt strong couphng, 



discussed in [|^, |T^, Here B\ B *, B" denote the relevant B meson poles, 
and /i^2 denotes the pseudoscalar meson decay constants. The coupling /3 
has been analyzed in and found to be close to zero and therefore will be 
neglected. 

The operator O4 has the same kind of decomposition as Oi, while the 
operator Oq can be rewritten as the product of density operators. The evalu- 
ation of the matrix elements like < K\q{l+'y5)b\B > and < 7r7r|g(l— 75)6|i? > 
can then be reduced to the evaluation of the matrix elements of the weak 
currents < 7r7r|g7^756|_B > and < Klq^y^blB > following the procedure de- 
scribed in detail in For the s(l + 75)^ scalar and pseudoscalar quark 
density operator we use the CHPT result and obtain 



s{l+^,)q = -S.BUl (18) 

where i3 is a real constant expressed in terms of quark and meson masses; 
e.g., to lowest order mj^o = B{ms + rrid) and U = exp{i2Il/ f) where 11 is a 
pseudoscalar meson matrix We use the value = 150 MeV for the s 
quark mass at the scale mass [H], giving B = 1.6 GeV. For the calculation 



of the density operator q{l — 75)6 we use the relations 

5756 = — 9„(g7"756), (19) 

and 

qb = —d^iqrb), (20) 

TUb 

where has been dropped since << rrib. 

The factorization assumption in the case of the Oq operator results in: 

<K'{p3)7:+{p,)7r-{p2)\Oe\B{pB)> = 
-2 < ^|s(l + 75)m|0 >< 7r7r|M(l -75)6|5 > . (21) 

u,d,s,c,b 
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Since the density operator in < TrK\s{l + 75)(i|0 > can be related to the 
matrix element of the corresponding current, and the form factors describing 
the matrix element of the current operator are dominated by light vector 
meson resonances, contributions from terms of the type < 7ri^|s(l + 75)c/|0 > 

< 711^(1 — 75)615 > is suppressed in the high energy region for the nonres- 
onant amplitude, and is neglected here for our purpose. There is also the 
contribution < i^7r7r|s(l + 75)u|0 > < 0|u(l — 75)6|-B > and the pole contri- 
bution of the type < KittiICsIK^pb) > (m| - m|-)~^ < K\s{l + 75)m|0 > 

< 01^(1 — 75)6|i? >. Here Cg is the chiral Lagrangian for the light pseu- 
doscalar mesons. These two terms give negligible contribution to the branch- 



ing ratio as expected from annihilation graphs in p3|. The contribution 
< K7ni\s{l — 'j5)b\B~ > < 0|s(l + 75)s|0 > does not contribute due to the 
momentum conservation. 

Using the s and t variables we derive: 

< K-{p3)n+{pi)7T-{p2)\Oe\B{pB) >nr = 



13 r /1/2 2 nr flh/ 2 2 \ 

1— /— — m^r — (mo — rrio — s)w 

rriB Is Is 



nr 



h 



s + 2t-m%-2m{- mDw"]^}. (22) 



The dominant contributions to the amplitude for i? K tt+tt decays are 
then 

G 

Mnr = ^[VubV*,ai < Knn\Oi\B >„^ -VtbVZiai < KmT\Oi\B >„^ + 



V2' 



ae < K7rn\0e\B >,,)]. (23) 



The decay width of the nonresonant B ^ vr n^K can be found with the 
help of 

KAB- - = ^^J\M„A^ As it. (24) 

In the case of B~ K~K^K~ one has to take into account the two identical 
mesons in the final state (factor 1/2 in the branching ratio) and to include 
in the amplitude the additional contribution obtained by the replacement 
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The lower and the upper bounds are given by Smin = {^i +^2)^, Smax = 
{rriB — Tn^Y, while for t they are given by 

tmin,max{s) = + TTlg - ^ [(s - m| + mg) (s + 7712 - m^) 

T \Hs,ml,ml)\^{s,ml,ml)], (25) 

where A (a, b, c) = + 1^ + (? — 2{ab + ac + ah). 

The Wolfenstein parametrization of the CKM matrix gives Vub = 
AX^ip - 7//), Vus = A, K« = 1 - AV2, Kfe = AX', Vts = -AX' and Vtb = 1. 
The values A = 0.82, A = 0.224, are well known, while for p and f] there are 



bounds given in |30|. This parametrization indicates Vub contains the phase, 
while the combination of V^*Vtb has none. We denote Vub = \Vub\e^'^ following 
[|]. In our numerical calculation we use the recent value /b = 175 MeV (see 
e.g. pll]). The strong coupling g is determined by the CLEO experimental 



results on the D* decay width [^. The value g = 0.57 reproduces this 



decay width. Heavy quark symmetry requires that this parameter be the 
same for B mesons. However, it has been pointed out that the off- mass-shell 



effects are very important in view of the higher B meson energy scale [33 



and therefore we will use the lower value oi g = 0.23 ^, This smaller 
value of g correctly reproduces the form factor for the B ^ n transition at 



zero momentum transfer, calculated in this model, |]21| in agreement with 
results obtained in other approaches (e.g. |^). We fix the parameters 
ai_2 to repropduce the value |y4f^*(0)| in |3^. These results are obtained 



by assuming a pole behaviour of form factors. We follow this procedure 
and by using the values of these form factors at Q^ax^ '^^ obtain af^* = 
0.16[GeV^/^] and af^* = 0.05[GeV^/% which after soft scaling described 
in [H, [1^ give for the B mesons af^* = af^ and a2^* /niD = 012'^ /niB 
resulting in the values ai = 0.16 02 = 0.15 [GeV^/^] which we will 

use here. These values give 2(0) for the B —y p transiton in agreement 
with the values used in [^. The branching ratio can be written as a sum of 
tree level contribution T (in which the operator Oi contributes), the penguin 
contributions P (operators O4 and Og) and there are two interference terms 
proportional to cos 7 (denoted by Ji) and sin 7 (denoted by 12): 

BR{B- K-Ti+Ti^)nr = T + P + Ji COS 7 + sin 7. (26) 

After performing the numerical integrations using the parameters described 
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above, we obtain T = 7.0 x 10"^ P = 7.5 x 10^^ h = -4.3 x 10"^ and 
/2 = -1.5 X 10"^ 

In the case of B- K-R+R- decay we find for T = 3.4 x 10"^ P = 
3.7 X 10-^ Ji = -2.1 X 10"^ and h = -7.4 x 10"^ Taking a larger value 



g ~ 0.4, as some of typical smaller value of g as mentioned in we find 
moderate increase of all values T, P, I12 by about factor 20%. Note that 
while the penguin contributions dominate in these decays, the size of 7 is 
very important for the magnitude of the branching ratio. 



The results obtained by Belle collaboration [|I^ are BR(S — * R vr+vr )exp 
= (55.6 ± 5.8 ± 7.7) x 10"^ and BR{B- R-R+R-)exp = (35.3 ± 3.7 ± 
4.3) X 10~^. Comparing our results with these data we conclude that larger 
values of g (e.g. g ~ 0.4) should be discounted in these cases; they give rates 
for nonresonant decays larger than the experimental results. 

Due to the strong phases in 04 and Og one can generate the CP violating 
asymmetry as 

^ \nB- ^ TT+TT-K-) - r(i?+ 7r+7r-ir+)| 

\r{B- ^n+7i-R-) + r{B+ ^n+n-R+)\' ^ ' 

which can be written in the form 

sin -yNi 

A = — — , (28) 

N2 + cos 7A^3 ^ ^ 

where 

iVi = -2G^\V:,\\V^t\\V*\\Vtb\ J {dPS)a, <0^> [\m{a^) < O4 > + 

Im(a6) < Oe >], (29) 



No 



G^\V:f\Vub\^ JidPS){ai\ < Oi > \^ + \V;f\Vtt\^[\ai\^ < O4 > + 

\a^\^<Oe >]}, (30) 

= -2G^\V:,\\Vub\\V*\\Vtb\ J {dPS)ai < Oi > [Re{a,) < O4 > + 

Re{ae) < >]. (31) 
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The numerical calculation gives for the B tt+tt K decay A^^i = 
-3.0 X 10-^ N2 = 16,4 X 10-5 and N3 = -8.6 x 10"^ For the ^ 
K+K-R-decay it becomes Ni = -1.5 x 10-^ N2 = 8.2 x 10"^ and N3 = 
—4.2 X 10"^. One might think that the nonresonant decay amplitude might 
have more complicated structure than what we found within this model and 
then the CP violating asymmetry will be changed. However, in the fac- 
torization model, which we use in the above calculation, the CP violating 
asymmetry in the nonresonant partial or integrated decay rates, is essentially 
model independent, since the matrix elements of O4 and Og are practically 
proportional to that of Oi, as can be seen from Eqs ( jl?]) and (p2[). Even in 
the presence of low-energy resonances contribution in the nonresonant ampli- 
tude, the CP asymmetry, for the total decay rates integrated over the whole 
phase space, or for the differential decay rates measured as a function of the 
two-pion or two-kaon invariant mass, can be computed in terms of the Wil- 
son coefficients 04 and ag, and are independent of the form factors. The CP 
violating asymmetry Eq.(^) depends on the strong phase generated by the 
absorptive part of the Wilson coefficients. 

On the other hand, the partial widths asymmetry at the charmonium res- 
onance discussed in |, |[ seems to be more reliable for a further constraint 
of the CP violating phase 7. 

In order to obtain the partial width CP asymmetry, one also needs to 
calculate the resonant decay amplitude B~ XcoK~ -n^n~K~ . This 
amplitude can be easily determined using the narrow width approximation, 
as in [0: 

M{B- ^ XcoK-) Mixco ^ vr+TT"). (32) 

In our numerical calculations we will use the recent Belle measurement of 
the branching ratio ^ BR(S± ^ Xco K^) = (6.0 ±1.1) x 10"^ The 
amplitude is then M.{B^ XcoK~) = 3.37 x 10~^ GeV. The Xco decay data 
then fix the decay amplitudes for \J^{XcO = 0.113 GeV and 

|A^(Xco ^ K-K+)\ = 0.126 GeV. 

The partial decay width Tp for B' MMK', M = 7r+, K+, which 
contains both the nonresonant and resonant contributions, is obtained then 
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by integration from Smin = ("^Xco " 21^^ J ^ g^^^ = (^^^^^ _^ 21^^0)2^ where 
m^^^ = 3.415 GeV and T^^^ = 0.0149t|]:[|[]l GeV is the width of the Xco- 

\ \ f ^max r imax (^) 

= ds dt \Mnr + -Mr-r. (33) 

Similarly, Tp, the partial decay width for 5+ ^ MMK+, M = K+ 
also contains both the nonresonant and resonant contributions. The CP 
nonconserving asymmetry is defined by 

A, = (34) 

In calculation of the Tp — Tp we derive, asumming that Vub = |Kib|e*'^ : 

r„ — = Sm7- ^ 



/Smax pirrLaxyS) 
ds / dt —/=\Vub\\V*g\ai < K'mi\Oi\B >nr 



(35) 



The denominator is given by 

^ ^ ~ (27r)332m| Js 

x{|A^„,p + \M{B- ^ Xco^-) 7^A^(XcO ^ vr-7r+)p}. 

(36) 



Then we can write 



__4isino_ 
^^-^2 + A3 cos 7' ^ ^ 



where Ai is determined by (|35D, while A2 contains the sum of the resonant 
decay amplitude as well as 7-independent part of (^), corresponding to N2 
with the integration over the Xco resonance region, Smin = {j^xco ~ 
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Smax = {f^xco + '^^Xc^o)^)- '^^^ P^^^ which is 7 dependent is given in A3 (cor- 
responding to A^3 with the integration over the Xco resonance region). Note 
that the term proportional to cos 7 arising from the interference with the res- 
onance in the denominator does not contribute. Namely, in the integration 
region of s in (|36|) , the real part of the resonance amplitude, being antisym- 
metric in s gives no contribution to the integrated decay rates over the Xco 
resonance. 

In Fig.l and Fig.2 we present the distribution dVp/ds in the region 
[smin, Smax]- It is clcar from these figures that the nonresonant contribution 
is approximately constant over this region. For the given set parameters, 
with g = 0.23, the CP violating asymmetries are 

MB^ K^n^n-) = ^'^'f^ , (38) 

73 — 1.2 cos 7 

and 

AJB^ ^ K^K+K-) = ^-^'^^^ (39) 

' 41-5.6COS7 ^ ^ 

Following recently given bound for the phase 7 = (58 ± 24)° |]l[] we find that 
the partial widths asymmetry for the B^ K^n^n^ decay is reaching 10% 
and 16% for the B^ K^K+K-. 

The uncertainties due to the errors in the remaining input parameters 
have not been included here, but we can roughly estimate that the error in 
the asymmetry can be as large as 40%. 

In summary, we can compare the previously considered cases B^ — > 
MMn^, M = n+,K+ [|, § with the present cases B^ MMK^, M = 
TT'^,K'^. Due to the different CKM matrix elements the penguin contribu- 
tions dominate in the nonresonant B"^ MMK^, M = 7r+, decay rate. 
Moreover, the phase 7 crucially influences their magnitudes. The appearence 
of the strong phases in the 04 and Og Wilson coefficients of the penguin op- 
erators can generate the CP violating asymmetry in the nonresonant decay 
modes B^ MMK^, M = 7r+, K+. 

The CP violating partial widths asymmetry is smaller in the cases B^ 
MMK^, M = Tr+,K+ than S± ^ MMn^, M = Tr+,K+. One might 
expect that, since now the CP violating partial widths asymmetry comes from 
the interference of the resonant amplitude and the tree-level nonresonant 
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amplitude, which is given by VubV*^ which is CKM suppressed, relative to 
KfeKd for the decays 5± MMvr^, M = 7t+,K+. One could also measure 
the differential decay rates and CP asymmetry as a function of the two-pion 
and two-kaon invariant mass squared as shown in Fig.l and Fig. 2 to obtain 
more informations on the CP asymmetry as well as possible determination 
of the nonresonant amplitudes at the Xco mass region. 

The determination of these partial width asymmetries provides useful 
guidance for experimental searches of CP violating effects and further con- 
straints on the phase 7. 
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Figure 1: Differential branching ratios for B — >■ tt+tt K vs. s. The curves 
(a), (b), (c) are rf(B)(NR)/ds, d{B)/ds + d{B)/ds, d{B) / ds - d{B) / ds against 
the two-pion invariant mass squared s, for the nonresonant, CP symmetric 
and CP antisymmetric differential branching ratios respectively. 
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Figure 2: Differential branching ratios for i? K^K K vs. s. The curves 
(a), (b), (c) are c?(B)(NR)/ds, d{B) / ds + d{B) / ds , d{B) / ds - d(B) / ds against 
the two-kaon invariant mass squared s, for the nonresonant, CP symmetric 
and CP antisymmetric differential branching ratios respectively. 
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